July, 2003 



Hadronic B Decays Involving Even Parity Charmed Mesons 

Hai-Yang Cheng 

Institute of Physics, Academia Sinica 
Taipei, Taiwan 115, Republic of China 

cn" 

O 

o 

CN ■ Abstract 



(N 



Hadronic B decays containing an even-parity charmed meson in the final state are stud- 
ied. Specifically we focus on the Cabibbo-allowed decays B -» D**ir(p), D**D^\ D**D^ and 
B s — ► D**ir(p), where D** denotes generically a p-wave charmed meson. The B — > D** transition 
form factors are studied in the improved version of the Isgur-Scora-Grinstein-Wise quark model. 
We apply heavy quark effective theory and chiral symmetry to study the strong decays of p-wave 
charmed mesons and determine the mag nitude of the D\ ,2 -E% 2 mixing angle (the superscript 
standing for the total angular momentum of the light quark). Except the decay to -Di(2427)°7r~ 
the predictions of B(B~ — > D**°ir~) agree with experiment. The sign of the D^ 2 — D^ 2 mix- 
ing angle is found to be positive in order to avoid a severe suppression on the production of 
-Di(2427)°7r~. The interference between color-allowed and color-suppressed tree amplitudes is ex- 
pected to be destructive in the decay B~ — ► Di(2427)°7r~. Hence, an observation of the ratio 
Z?i(2427) 7r~/Di(2427) + 7r~ can be used to test the relative signs of various form factors as implied 
by heavy quark symmetry. Although the predicted B~ — > L>i(2420)°p~ at the level of 3 x 10~ 3 
exceeds the present upper limit, it leads to the ratio -Di(2420)/9~/Di(2420)7r~ 2.6 as expected 
from the factorization approach and from the ratio fp/fn ~ 1.6. Therefore, it is crucial to have 
a measurement of this mode to test the factorization hypothesis. For B — > D** D decays, it is 
expected that D* s0 D > D S \D as the decay constants of the multiplet (D* Q ,D s x) become identical 
in the heavy quark limit. The experimental observation of less abundant production of D* Q D than 
D S \D is thus a surprise. What is the cause for the discrepancy between theory and experiment 
remains unclear. Meanwhile, it is also important to measure the B decay into D s i(2536)-D(*) 
to see if it is suppressed relative to Z) s i(2460)-D(*) to test the heavy quark symmetry relation 
/d s i(2536) ^ /d s i(2460)- Under the factorization hypothesis, the production of D* s2 D is prohibited 
as the tensor meson cannot be produced from the V — A current. Nevertheless, it can be induced via 
final-state interactions or nonfactorizable contributions and hence an observation of B — * D s2 D^*' 
could imply the importance of final-state rescattering effects. 
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I. INTRODUCTION 



The interest in even-parity charmed mesons is revived by a recent discovery of a new 
narrow resonance by BaBar [1]. This state which can be identified with J p = + is lighter 
than most theoretical predictions for a + cs state [2]. Moreover, a renewed lattice calculation 
[3] yields a larger mass than what is observed. This unexpected and surprising disparity 
between theory and experiment has sparked a flurry of many theory papers. For example, it 
has been advocated that this new state is a four-quark bound state [4,5],* or a DK molecular 
[9] or even a Dir atom [10]. On the contrary, it has been put forward that, based on heavy 
quark effective theory and chiral perturbation theory, the newly observed £) s (2317) is the 
+ cs state and that there is a 1 + chiral partner with the same mass splitting with respect 
to the Y~ state as that between the + and 0~ states [11,12]. The existence of a new narrow 
resonance with a mass near 2.46 GeV which can be identified with 1 + state was first hinted 
by BaBar and has been observed and established by CLEO [13] and Belle [14]. 

Although the -D* (2317) and -D s i(2460) states were discovered in charm fragmentation 
of e + e~ — > cc, it will be much more difficult to measure the counterpart of _D* (2317) and 
L> s i(2460) in the non-strange charm sector, namely D$ and Di, owing to their large widths. 
Indeed, the broad Dq and D\ resonances were explored by Belle [14] in charged B to D + 7t~tt~ 
and D* + 7t~tt~ decays (see Table I). The study of even-parity charmed meson production in 
B decays, which is the main object of this paper, also provides an opportunity to test heavy 
quark effective theory. 

This work is organized as follows. The masses and widths of p-wave charmed mesons are 
summarized in Sec. II. In order to determine the mixing angle of the axial-vector charmed 
mesons, we apply heavy quark effective theory and chiral symmetry to study their strong 
decays. The decay constants of p-wave charmed mesons and B — > D** form factors are 



*The low- lying non-charm scalar mesons in the conventional qq' states are predicted by the quark 
potential model to lie in between 1 and 2 GeV, corresponding to the nonet states /o(1370), do (1450), 
#5(1430) and / (1500)// (1710). The light scalar nonet formed by cr(600), k(800), /o(980) and 
ao (980) can be identified primarily as four-quark states [6]. It has been argued [7] that a strong 
attraction between (qq)s* and (qq)s [6,8], where 3* and 3 here refer to color, and the absence of the 
orbital angular momentum barrier in the s-wave 4-quark state, may explain why the scalar nonet 
formed by four-quark bound states is lighter than the conventional qq nonet. By the same token, 
it is likely that a scalar cnns 4-quark state, where n = u, d, will be lighter than the + p-wave cs 
state, where a typical potential model prediction gives 2487 MeV [2]. It has been suggested in [4] to 
search for exotic 4-quark cqqq charmed meson production in B decays. Particularly noteworthy are 
resonances in the doubly charged Dfir + (D + K + ), and wrong pairing D + K~ channels. However, 
contrary to the case of scalar resonances, the 1 + D s (2460) state is unlikely a four-quark state as it 
is heavier than the axial-vector meson formed by cs. A non-observation of a heavier and broad + 
cs state will not support the four-quark interpretation of D s (2317). 
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studied in Sec. Ill within the Isgur-Scora-Grinstein-Wise quark model. The production of 
p-wave charmed mesons in B decays is studied in detail in Sec. IV. Conclusions are presented 
in Sec. V. 

TABLE I. The masses and decay widths of even-parity charmed mesons. We follow the naming 
scheme of Particle Data Group [15] to add a superscript "*" to the states if the spin-parity is in 
the "normal" sense, J p = + , 1~ , 2 + , • • •. The four p-wave charmed meson states are thus denoted 
by Dq,Di,D[ and D\. In the heavy quark limit, D\ has j = 1/2 and D[ has j = 3/2 with j being 
the total angular momentum of the light degrees of freedom. 



State 


Mass (MeV) 


Width (MeV) 


Reference 


£>5(2308)° 


2308 ± 17 ± 15 ±28 


276 ± 21 ± 18 ± 60 


[14] 


£>i(2427) 


2427 ± 26 ± 20 ± 15 


384 +1 °J ± 24 ± 70 


[14] 


£>i(2420)° 


2422.2 ± 1.8 


i8.9l|;| 


[15] 




2421.4 ± 1.5 ±0.4 ±0.8 


23.7 ±2.7 ±0.2 ±4.0 


[14] 


Z? / 1 (2420) ± 


2427 ± 5 


28 ±8 


[15] 


1^(2460)° 


2458.9 ± 2.0 


23 ±5 


[15] 




2461.6 ±2.1 ±0.5 ±3.3 


45.6 ± 4.4 ± 6.5 ± 1.6 


[14] 


^(2460) ± 


2459 ± 4 


25l 8 7 


[15] 


£> s *o(2317) 


2317.3 ±0.4 


< 7 


[1,13,14] 


£ s i(2460) 


2458.9 ± 1.2 


< 7 


[13,14] 


^i(2536) 


2535.35 ±0.34 ±0.5 


< 2.3 


[15] 


L>* 2 (2573) 


2572.4 ± 1.5 


15±t 


[15] 



II. MASS SPECTRUM AND DECAY WIDTH 

In the quark model, the even-parity mesons are conventionally classified according to the 
quantum numbers J, L, S: the scalar and tensor mesons correspond to 2S+1 Lj = 3 P and 
3 P 2 , respectively, and there exit two different axial-vector meson states, namely, 1 P\ and 
3 Pi which can undergo mixing if the two constituent quarks do not have the same masses. 
For heavy mesons, the heavy quark spin Sq decouples from the other degrees of freedom in 
the heavy quark limit, so that Sq and the total angular momentum of the light quark j are 
separately good quantum numbers. The total angular momentum J of the meson is given 
by J — j + Sq with S = s + Sq being the total spin angular momentum. Consequently, it is 
more natural to use = P^ 2 , P^ 2 ■, P^ 2 and P^ 2 to classify the first excited heavy meson 
states where L here is the orbital angular momentum of the light quark. It is obvious that 
the first and last of these states are 3 P 2 and 3 P , while [17] 
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In the heavy quark limit, the physical eigenstates with J p = 1 + are Pj and P 1 / rather 
than 3 Px and 1 P 1 . 

The masses and decay widths of even-parity (or p-wave)^ charmed mesons Dj and D* sJ 
are summarized in Table I. We shall use 1 + and 1 /+ or D\ and D[ to distinguish between 
two different physical axial-vector charmed meson states. As we shall see below, the physical 
f + state is primarily P^ 2 , while f /+ is predominately P^ 2 . A similar broad Di state not 
listed in Table I was reported by CLEO [f 8] with M = 2461±i MeV and T = 2901^ MeV. 
For the known narrow resonances, the Belle's measurement of the D 2 ° width (see Table I) 
is substantially higher than the current world average of 23 ± 5 MeV [15]. 

In the heavy quark limit, the states within the chiral doublets (0 + , 1 + ) with j = 1/2 and 
(1 /+ , 2 + ) with j = 3/2 are degenerate. After spontaneous chiral symmetry breaking, + states 
acquire masses while 0~ states become massless Goldstone bosons. As shown in [11], the fine 
splitting between + and CT is proportional to the constituent quark mass. The hyperfine 
mass splittings of the four p-wave charmed meson states arise from spin-orbit interactions 
and tensor-force interaction [see Eq. (2.21) below], while the spin-spin interaction is solely 
responsible for the hyperfine splitting within the multiplet (CT, 

From Table I and the given masses of pseudoscalar and vector charmed mesons in the 
PDG [15], it is found empirically that the hyperfine splitting within the chiral multiplets 
(0 + , 1 + ), (l /+ ,2 + ) and (0~, 1~) are independent of the flavor of the light quark: 

m(D* 2 ) - m(D[) « m(D* 2 ) - m(D' sl ) « 37 MeV, 
m(D 9l ) - m(D* s0 ) « m(D* s ) - m{D s ) = 144 MeV, 

m{D 1 ) - m(D*) w m(D*) - m(D) = 143 MeV. (2.2) 
However, the fine splittings 

m(D* s0 ) - m(D s ) w m(D sl ) - m(D* s ) sa 350 MeV, 

m(D*) - m{D) w m{D l ) - m(D*) w 430 MeV, (2.3) 

depend on the light quark flavor. * Since the fine splitting between + and 0" or 1 + and 
1~ should be heavy-flavor independent in heavy quark limit, the experimental result (2.3) 



tlf the even-parity mesons are the bound states of four quarks, they are in an orbital s-wave. In 
this case, one uses J p = + rather 3 P to denote scalar mesons, for example. 

■I-Likewise, considering the spin-averaged masses of the doublets (0 + , 1 + ) and (l /+ ,2 + ) 

moi(D) = \m(Dl) + |m(Di), m 12 (D) = ^m(D[) + ^m(D* 2 ), (2.4) 

the hyperfine mass splittings 

m 12 (D) -moi(-D) w 48MeV, m 12 (D s ) - m 01 (D s ) w 132MeV, (2.5) 
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implies that the fine splitting is light quark mass dependent. Indeed, if the first line rather 
than the second line of Eq. (2.3) is employed as an input for the fine splittings of non-strange 
charmed mesons one will predict [11] 

M(Dq ± ) = 2217 MeV, M(D* °) = 2212 MeV, 

M{Df) = 2358 MeV, M(£>°) = 2355 MeV, (2.6) 

which are evidently smaller than what are measured by Belle [14]. 

It is interesting to note that the mass difference between strange and non-strange charmed 
mesons is of order 100 ~ 110 MeV for 0~, 1 _ , 1 /+ and 2 + as expected from the quark model. 
As a consequence, the experimental fact that m(D* s0 ) iti(Dq) and m(D sl ) ~ m(Di) is very 
surprising. 

In the heavy quark limit, the physical mass eigenstates D\ and D[ can be identified 

1/2 3/2 

with P 1 and P 1 / , respectively. However, beyond the heavy quark limit, there is a mixing 
between P^ 2 and P^ 2 , denoted by D^ 2 and D^ 2 respectively, 

D 1 (2427) = D\ /2 cos 9 + D\ /2 sin 6, 

L>i(2420) = -D\ /2 sin 9 + D\ 12 cos 9. (2.7) 
Likewise for strange axial-vector charmed mesons 

D sl (2460) = D][ 2 cos9 s + D 3 s ( 2 sm9 s , 

D' sl (2536) = —D X J 2 sin 9 S + Dl[ 2 cos9 s . (2.8) 
Since D^ 2 is much broader than D^ 2 as we shall see shortly, the decay width of D^(2420) 

1/2 3/2 

is sensitive to the mixing angle 9. Our task is to determine the D ± —D^ mixing angle 
from the measured widths. In contrast, the present upper limits on the widths of _D s i(2460) 
and D' sl (2536) do not allow us to get any constraints on the mixing angle 9 S . Hence, we will 
turn to the quark potential model to extract 6 S as will be shown below. 

It is suitable and convenient to study the strong decays of heavy mesons within the 
framework of heavy quark effective theory in which heavy quark symmetry and chiral sym- 
metry are combined [21]. It is straightforward to generalize the formalism to heavy mesons 
in p-wave excited states [22]. The decay Dq undergoes a s-wave hadronic decay to Dn, while 

1/2 

Di can decay into D* by s-wave and ci-wave pion emissions but only the former is allowed 
in the heavy quark limit m c — > oo: 

T ( D^D« ) = A . D ^, r^ir,)^^*-, (2.9) 



also depend on the light quark flavor. Based on a quark-meson model, the spin-weighted masses 
mo-iip) = 2165 ± 50 MeV [19] and m 01 p s ) = 2411 ± 25 MeV [20] were predicted, while experi- 
mentally they are very similar (see Table I). 
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where p c is the cm. momentum of the final-state particles in the B rest frame. The tensor 
meson D\ decays into D* or D via ci-wave pion emission. In the heavy quark limit where 
the total angular momentum j of the light quark is conserved, D^ 2 — > Dtt is prohibited by 
heavy quark spin symmetry. The explicit expressions for the decay rates are [22], 

F(D 1 ^ D Ti) = — — , 

6tt m D 3/2 A 2 / 2 



IOtt m D2 A 2 /2' 

1 mo /i' 2 
157T m fl2 A 2 / 2 : 



r(D* - Dtt) = — — % (2.10) 



where A x is a chiral symmetry breaking scale, f n = 132 MeV and h' is a heavy flavor 
independent coupling constant. The pi dependence of the decay rate indicates the ci-wave 
nature of pion emission. From Eq. (2.10) we obtain 



T(Df -> D+iT-) 2 m D ( p c (D* -> £> 



7T 
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r(D|° -> D*+7r-) 3 m D * \p c {D* 2 ^ D*7i) 



= 2.3, (2.11) 



in excellent agreement with the measured value of 2.3 ± 0.6 [15]. 

Since the ci-wave decay is severely phase-space suppressed, it is evident that Dq and D\ 
are very broad, of order 250 MeV in their widths, whereas D[ and D\ are narrow with widths 
of order 20 MeV. 

The strong couplings appearing in Eq. (2.9) are given by 

m D - m D h 



9D*Dn = y/m Do m D °- — , 

m D Jit 



u i u n V u i m D i/2 U 



2 2 

m n i/2 — rn Dt . 7 

(2.12) 



with h being another heavy-flavor independent coupling constant in the effective Lagrangian 
[22]. It can be extracted from the measured width of -Dp (2308) ( see Table I) to be 

h = 0.65 ± 0.12 . (2.13) 

From the averaged width 29.4 ± 4.2 MeV measured for D^ we obtain 

h! 

— = 0.67 ±0.05 GeV" 1 . (2.14) 

Substituting the couplings h and h! into Eqs. (2.9) and (2.10) leads to 

T(Dl /2 -> D*tt) = 10.5 ± 1.5 MeV, Y{D\^ 2 — > D*tt) = 248 ± 92 MeV, (2.15) 
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where we have assumed that D x ' has a mass close to £)i(2420) and m(D{' ) « m(L>i(2427)). 

3/2 1/2 

Therefore, is much narrower than D ± ' owing to the phase-space suppression for ci-wave. 
However, the physical state D' 1 (2A20) can receive a s-wave contribution as there is a mixing 

1 /2 3 /2 

between D± and Z) 1 beyond the heavy quark limit. The observed narrowness of D 1 (2420) 
indicates that the mixing angle should be small; that is, £^(2420) should be dominated by 
D^ 2 while -Di(2427) is primarily the D^ 2 state. Using 

r[/J / 1 (2420)] = V{D\ /2 -> D*7i) cos 2 9 + r(D 1 1/2 -> D*n) sin 2 (2.16) 

and the averaged width 20.9 ± 2.1 MeV for D' 1 (2420)°, it is found § 

e = ±(12. i+6-6)°. Q.n) 

We shall see in Sec. IV that a positive mixing angle is preferred by a study of D 1 (2427)°ir~ 
production in B decays. 

The scalar resonance -D* (2317) is below the threshold of DK and its only allowed strong 
decay -D* (2317) — > D s n° is isospin violating. Therefore, it is extremely narrow with a width 
of order 10 keV [4,11,23]. As for D s i(2460), it is below D*K threshold and its decay to 
DK is forbidden by parity and angular momentum conservation. Hence, the allowed strong 
decays are D*7r°, D* 7r°, D s 7ttt and D s 7T7ttt. The isospin- violating decay -D s i(2460) — > D*7r° 
has a rate similar to -D* (2317) — > D s ir°. At a first sight, it is tempted to argue that D s im 
could dominate over D*tt° as the former can proceed without violating isospin symmetry. 
However, a detailed analysis shows this is not the case. The decay D sl — > D s iiir arises from 
the weak transitions -D s i(2460) — > D s {a(600), / (980), • • •} followed by the strong decays 
{a"(600), /o(980), • • •} — > 7T7T. Consider the dominant contributions from the intermediate 
states <t(600), / (980) and their mixing 

/o = ss cos + nn sin 0, a = — ss sin + nn cos 0, (2-18) 

with nn = {uu + dd) /y/2. The u — fo mixing angle can be inferred from various processes, 
see [24] for a summary. In general, the mixing angle is small so that /o(980) has a large 
ss component while a is primarily nn. The /o production is favored by the weak decay of 
D s i into D s , but its contribution to tttt is suppressed by the large off-shellness of /o(980), 
recalling that the mass difference between D s i(2460) and D s is only 494 MeV. In contrast, 
(j(600) is favored by phase space consideration and yet its contribution is suppressed by 
the small a — f Q mixing angle. As a net result, although the strong decay into D s nii is 
isospin conserving, its OZI suppression is more severe than the isospin one for D* s ir°. This 
is confirmed by the recent measurement of CLEO [13] 



n 1/2 3/2 

mixing angle was just reported to be 9 = 0.10 ± 0.03 ± 0.02 ± 0.02 radians = 
(5.7 ± 2.4)° by Belle through a detailed B — > D*irir analysis [14]. This is consistent with the result 
ofEq. (2.17). 
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g(^i(2460) -» D 
B(D sl (2AQ0) ^ D* s ^) <0 -° 8 - (2 ' 19) 

As for the electromagnetic decays of £> sl (2460), CLEO and Belle found 



g(D al (2460) -> £> a7 ) 
B(D sl (2460) -> D*7T°) 



< 0.49 CLEO [13], 

0.38 ±0.11 ±0.04 Belle [14], (2.20) 

0.63 ±0.15 ±0.15 Belle [16]. 



There are two different measurements of the radiative mode by Belle: the measured value 
0.38 ± 0.11 ± 0.04 is determined from B — > DD sl decays, while the other value comes from 
charm fragmentation of e + e~ — > cc. These two measurements are consistent with each 
other, though the central values differ significantly. Hence, just as its + partner D* (2317), 
-D s i(2460) is also extremely narrow. A theoretical estimation yields 38.2 keV for its width 
[11]. 

Eq. (2.14) leads to T(D* 2 ) = 12.6 MeV, in agreement with experiment (see Table I). 

3 /2 1 /2 

Since T(D sl ) = 280 keV followed from Eqs. (2.10) and (2.14) and D sl is very narrow as 
its mass is close to -D sl (2460) which is below D*K threshold, the decay width of D' sl (2536) 
is thus at most of order 0.3 MeV and is consistent with the experimental limit 2.3 MeV [15]. 
In short, while D[ and D\ are rather narrow, Dq and D\ are quite broad as they are allowed 
to have s-wave hadronic decays. In sharp contrast, D* and D s \ are even much narrower 
than D' sl as their allowed strong decays are isospin violating. 

Since the width of D' sl (2536) has not been measured, we will appeal to the quark potential 
model to estimate the Dl{ 2 — D 3 s { 2 mixing angle [see Eq. (2.8)]. It is known that spin-orbital 
and tensor force interactions are responsible for the mass splitting of the four p-wave charmed 
mesons. In the quark potential model the relevant mass operator has the form [25] 

M = \£ ■ S! + 4t£ ■ s 2 + rS 12 , (2.21) 

with si and s 2 refer to the spin of the light and heavy quarks, respectively, and 

1 



A 



2m? q 



r \ m c / r 
V 



Am q m c r 



(2.22) 



where V(r) is the zero-component of a vector potential, S(r) is a scalar potential responsible 
for confinement, and Su is the tensor force operator 

= 3<r i ■ r cr 2 ■ r — <t x ■ er 2 . (2.23) 

Note that the assumption of a Coulomb-like potential for V(r) has been made in deriving 
Eq. (2.21) [25]. Under this hypothesis, the mass splitting is governed by the two parameters 
A and r. 
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Following Calm and Jackson [25], the masses of the J = 2 and J = states read 

A 8 

M 2 = - + -r + c, M = -A - 8r + c, (2.24) 
2 5 

while the masses of the two J = 1 states obtained by diagonalizing the matrix in the 
\J,j, to) = |l,3/2, to) and |1, 1/2, to) bases are (up to a common mass c) 



A _ 8 _2a/| t . ~ 
2 3 3 



(2.25) 



It is clear that the mixing vanishes in the heavy quark limit m c — > oo. However, 1/to c 
corrections will allow charm quark spin to flip and mix D^ 2 and D^/ 2 . The two eigenmasses 
for J = 1 are then 



M 1± = -^ + c±J^ + i(A-4r)2. (2.26) 



From Eqs. (2.8) and (2.25) we arrive at 



= S i n -i I i?+ ) with i?+ = ^ * T Ml+ . (2.27) 



3 ' 



** 



The parameters A and r are obtained by a global fit to the charm spectroscopy. For D 
mesons, it is found 

r^l2MeV, A xs 104 MeV, 9 S w 7°. (2.28) 

As pointed out in [25], a positive spin-orbit energy A implies a less important scalar potential 
S. On the contrary, the existing potential model calculation such as the one by Di Pierro 
and Eichten [2] yields A < [25] or a very strong confining potential S. This will cause a 
reversed splitting, namely the j = 1/2 states lying above the j = 3/2 states. However, we 
will not address this issue here. For D** mesons, the parameters A and r fall into some large 
regions because of large uncertainties associated with the measured masses of Dq and D\. 
Hence, the magnitude and even the sign of the D\ l2 -Dl /2 mixing angle 6 at present cannot 
be fixed within this approach. Instead, we have used heavy quark effective theory together 
with the measured widths to extract \9\. As will be seen below, the sign of 6 can be inferred 
from a study of the p-wave charmed meson production in B decays. 

III. DECAY CONSTANTS AND FORM FACTORS 
A. Decay constants 

The decay constants of scalar and pseudoscalar mesons are defined by 
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(0\A,\P(q))=if P q„, 



(0\V,\S(q)) = f s q,. 



(3.1) 



It is known that the decay constants of non-charm light scalar mesons are smaller than that 
of pseudoscalar mesons as they vanish in the SU(3) limit. For the neutral scalars er(600), 
/o(980) and Oq(980), the decay constant must be zero owing to charge conjugation invariance 
or conservation of vector current: 

U = ffo = /a° = 0. (3.2) 

Applying the equation of motion, it is easily seen that the decay constant of Kq + (clq) is pro- 
portional to the mass difference between the constituent s (d) and u quarks. Consequently, 
the decay constant of the charged ao(980) is very small, while the one for 1^(1430) is less 
suppressed. A calculation based on the finite-energy sum rules [26] yields 

f a ± = 1.1 MeV, f K . = 42 MeV. (3.3) 

Contrary to the non-charm scalar resonances, the decay constant of the scalar charmed 
meson is not expected to be suppressed because of charm and light quark mass imbalance. 
Applying the equation of motion again leads to 

m K*fK* = i{m s ~ m u )(K^\su\0), m D*fD* = i(m c ~ m u ) (Dq\cu\0) . (3.4) 

For a crude estimate, we assume (Dq\cu\0) (Kq\su\0) and obtain 

fo* ~ 160 MeV. (3.5) 

This is comparable to f D 200 MeV, the decay constant of the pseudoscalar D meson. 
The decay constants of the axial-vector charmed mesons are defined by 

(Q\ A n\ D i 2 (l,£)) = f D v* m D w £ i» (0|-4 M |-Di /2 (?,£)) = / D 3/2m D 3 /2 e M . (3.6) 
It has been shown that in the heavy quark limit [27-29] 

/ D l/2 = fD*, f D 3/2 =0. (3.7) 

Since the decay constant of D\ vanishes irrespective of heavy quark symmetry (see below), 
the charmed mesons within the multiplet (0 + , 1 + ) or (l /+ ,2 + ) thus have the same decay 
constant. This is opposite to the case of light j9-wave mesons where the decay constant of 
1 Pi meson vanishes in the SU(3) limit [30] based on the argument that for non-charm axial 
vector mesons, the 3 Pi and l P\ states transfer under charge conjunction as 

M b a ( 3 P 1 ) - M«( 3 Pi), M b a CPi) - -M£( 1 P 1 ), (a = 1, 2, 3), (3.8) 

where the axial- vector mesons are represented by a 3 x 3 matrix. Since the weak axial- vector 
current transfers as (A^) b a — > {A^)l under charge conjugation, it is clear that the decay 
constant of the l P\ meson vanishes in the SU(3) limit [30]. 
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The polarization tensor of a tensor meson satisfies the relations 

e^ = e vlx , e% = Q, p^ v = p v e^ = 0. (3.9) 

Therefore, 

(0|(y - A)^\D*(e,p)) = ae^tf + be\ Vlx = 0. (3.10) 

The above relation in general follows from Lorentz covariance and parity consideration. 
Hence the decay constant of the tensor meson vanishes; that is, the tensor meson D\ cannot 
be produced from the V — A current. 

Beyond the heavy quark limit, the relations (3.7) receive large l/m c corrections which 
have been estimated in [29] using the relativistic quark model. In the present paper we shall 
use f p = 216 MeV and (in units of MeV) 

f D = 200, f Ds = 230, f D . = 230, 

f Do = 160, / fl i/ S = 120, / D 3/ 2 = 40, 

/d s0 = 60, / d i/ 2 = 145, / D 3/ 2 =50. (3.11) 

si si 

The decay constants for D** are obtained essentially by fitting to experiment to be discussed 
in Sec. IV. B. Note that the measurements of B — > £)(*)£)(*) [15,31] indicate that the decay 
constants of D* and D s are similar. 



B. Form factors 

Form factors for B — > M transitions with M being a parity-odd meson are given by [32] 
(P(p)\V»\B(p B )) = (ips+p), -^-^^ Frtf) + ^^S^V), 
(V(p,e)\V,\B(p B )) = \ 6„ ua pe*»p B pPV(q 2 ), 

(V(p,e)\A,\B(p B )) = iUm P + my^Atf) - e *' PB (p B + p),A 2 (q 2 ) 

p mp + m v 



-2m v —^q,[A 3 (q 2 )-A (q 2 )]}, (3.12) 



where q = p B - p, F 1 (0) = F (0), A 3 (0) = A (0), and 



2 m P + m v 2 mp-m v 2 
A 3 (g ) = A(<Z ) ^— A 2 (, ). (3.13) 

For i? — > P and I? — > V form factors, we will use the Melikhov-Stech (MS) model [33] based 
on the constituent quark picture. Other form factor models give similar results. 

The general expressions for B — > D** transitions (D** being a p-wave charmed meson) 
are given by [34] 
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(D*(p)\A„\B(p B )) = % u + {q 2 ){p B +p) IM + u_(q 2 )(p B 



{D\ /2 {p,s)\V,\B{p B )) 



'V 



2 (g 2 )< + cf • Ps )(p s + p)„ + ci' V)(e* • Ps )(p s - p) 



1/2, 



(£) 1 1/2 (p,e)|A M |S(p B )) = -qi/2(q 2 Ku P ae* u (pB + pY(Pb -pf, 



(Dl /2 (p,e)\V„\B(p B ))=i e V 2(q 2 )e; + cr(q 2 )(e* .p B )(p B +p), + cT(q 2 )(e* .p B )(p B -p), 



3/2 



3/2, 



(D? /2 (p,e)|^|S(p B )) = -g 3/2 (g 2 )e^ p(T ^(pB +pY(pb -pY, 
(D*(pe)\V,\B(p B )) = h(q 2 Y, upa e*» a (p B ) a (p B +pY(pB-pY, 
(D*(p,e)\A,\B(p B )) = -i[k(q 2 )e; u p B + b + (q 2 )e^p B p B (p B +p), 

+ b„{q 2 y ap p B p B {p B -p)» . 



(3.14) 



In order to know the sign of various form factors appearing in Eq. (3.14), it is instructive to 
check the heavy quark limit behavior of B — > D** transitions which has the form [17] 



(D*(v')\A,\B(v 
(D{ /2 (v',e)\V,\B(v 

(Dl /2 (v',e)\A^\B(v 



{D\ /2 {v\e)\V,\B{v 



(Dl /2 (v',e)\A„\B(v 

(d;(v',e)\v,\b(v 

(D*(v',e)\A^\B(v 



y/m B m Do 2r 1/2 (u)i(v' - v)^ 
m B m D i/2 2r l/2 (uj)i 



u-1)eI-{e*.v)v^ 

v l a ft 



^m B m D i /2 2r 1/2 (uj)(-)e^ a pe* v v' v p , 
J-m B m D 3/2 r 3/2 (w)i{(l - uj 2 '^ - (e* ■ v)[3v^ - (u - 2)^]}, 



-m B m D s/ 2 t 3/2 (uj)(uj + l)e^ a/3 £*V Q f /3 , 

^3m B m D2 r ? ,/ 2 (uj)e flval3 e* ul v 1 v ia v (3 , 
yj3m B m D2 T 3 /2(u})(-i){(u + l)s*X - £*^V^}, 



(3.15) 



where ou = v ■ v' and there are two independent functions Ti/ 2 (ou) and t 3 i 2 {u) first introduced 
in [17]. It is easily seen that the matrix elements of weak currents vanish at the zero recoil 
point uj = 1 owing to the orthogonality of the wave functions of B and D**. The universal 
functions Ti/ 2 (cu) and t 3 / 2 {uj) are conventionally parametrized as 



T 1 {U)=T 1 {1)[1-P 2 {U-1)} 



(3.16) 



for i = 1/2 and 3/2. The slope parameter p 2 can be related to T!/ 2 (l) and t 3 / 2 (1) via 
the Bjorken sum rule. Comparing Eq. (3.14) with Eq. (3.15) we see that heavy quark 
symmetry requires that the form factors u + , £i/ 2 , q±/ 2 , c]/ 2 , h, k and 6_ be positive, while 
■u_, £3/2, <?3/ 2 , C+ 2 , c+ 2 , ci 2 and b + be negative. Heavy quark symmetry also demands the 
relations c 1 / 2 + ci 2 = and b + + 6_ = 0. It is easily seen that these heavy quark symmetry 
requirements are satisfied in realistic model calculations shown below. 

In the present paper, we shall use the improved version, the so-called ISGW2 model [35], 
of the non-relativistic quark model by Isgur-Scora-Grinstein-Wise (ISGW) [34] to compute 
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the B — > D** transition form factors.** In general, the form factors evaluated in the original 
version of the ISGW model are reliable only at q 2 = q 2 ^, the maximum momentum transfer. 
The reason is that the form-factor q 2 dependence in the ISGW model is proportional to 
exp[— {*q 2 m — q 2 )} and hence the form factor decreases exponentially as a function of (q 2 ^ — q 2 ). 
This has been improved in the ISGW2 model in which the form factor has a more realistic 
behavior at large (g^ — q 2 ) which is expressed in terms of a certain polynomial term. In 
addition to the form-factor momentum dependence, the ISGW2 model incorporates a num- 
ber of improvements, such as the constraints imposed by heavy quark symmetry, hyperfine 
distortions of wave functions, etc.,- • • [35]. The results of the ISGW2 model predictions for 
various form factors are shown in Tables II-IV. Evidently, the signs of various calculated 
form factors are consistent with what are expected from heavy quark symmetry. 

TABLE II. The form factors at various q 2 for B —>■ Dq and B s — ► D* transitions calculated in 
the ISGW2 model. 



Transition u + (m 2 ) u^(m 2 ) u + {m 2 ) u^{m 2 p ) u + (m 2 Dg ) u_{m 2 Ds ) 

B -► D* 0.175 -0.462 0.178 -0.471 0.198 -0.524 

B s -► D* s0 0.196 -0.515 0.200 -0.527 0.230 -0.605 



TABLE III. The form factors at q 2 = ml for B — > D 1 ' and B — > £>/ transitions calculated 
in the ISGW2 model, where £ 1 / 2 and are in units of GeV and all others are in units of GeV -1 . 



Transition 


Ql/2 


k/2 


c 1/2 


1/2 

C_ 


Q3/2 


h/2 


3/2 


3/2 

C_ 


B - D\ 12 


0.057 


0.54 


-0.064 


0.068 










B -► d\ /2 










-0.057 


-1.15 


-0.043 


-0.018 


Bs - DT 


0.063 


0.66 


-0.072 


0.078 










r — v n 3//2 










-0.063 


-1.31 


-0.048 


-0.023 



In realistic calculations of decay amplitudes it is convenient to employ the dimensionless 
form factors defined by [32] 



(D*(p Do )\A,\B(p B )) 



(Pb+Pd )h ~ 



m% 



m 



D 



Q,)F 1 BD °(<1 2 ) + 



m. 



m 



D 



**Note that in the original version of the ISGW model [34], the form factors for B to axial- vector 
charmed meson transition are evaluated for Z?i( 1 Pi) and Z?i( 3 Pi). As a result, one has to apply 
Eq. (2.1) to obtain B — ► D\^ 2 and B — > D^ 2 form factors. 
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TABLE IV. The B D* 2 and B s -► L>* 2 form factors at g 2 = m 2 calculated in the ISGW2 
model, where k is dimensionless and h, b+, b- are in units of GeV~ 2 . 



Transition h k b + 6_ 

B -> D* 2 0.011 0.60 -0.010 0.010 

B s -► D* 2 0.013 0.70 -0.011 0.012 



TABLE V. The form factors F and F\ at various q 2 for B ^ Dq and i? s — > D* transitions 
calculated in the ISGW2 model. 



Transition 






^iK) 




Fi(m 2 D J 


F (m 2 D J 


B^D* 


0.175 


0.175 


0.178 


0.166 


0.198 


0.108 


Bs - D* s0 


0.196 


0.196 


0.200 


0.187 


0.230 


0.130 



(D 1 (p Dl ,s)\V,\B(p B )) = i\ (m B + mn^Vf^iq 2 ) - 8 ' PB (p B + p Dl ),V 2 BD ^q 2 ) 

- 2m Dl £ -^{p B -p Dl ), [V^tf) - Vr^q 2 )} }, 

(D 1 (p Dl ,e)\A fl \B(p B )) = 1 e^e*»p B p° Dl A BD i(q 2 ), (3.17) 

m B + m Dl 

with 

Vf^(q 2 ) = mB + m ^ Vl ^(q 2 ) _ ^ -m gl Bg 2 
2m Dl 2m Dl 

and Kj BDl (0) = V^ Dl (0). The form factors relevant for B -> Z^P decays are P BD ° and P BP 
Note that only the form factor V BP or Ff Dl will contribute to the factorizable amplitude of 
B — > .DiP as one can check the matrix elements q^(D 1 (p B>1 ,e)\V fl \B(p B )) and e tJ, {P\V fl \B). 
The ISGW2 model predictions for the form factors F 0i , Vo,i,2 an d ^4 are summarized in 
Tables V-VI. It is evident that the form factor Po,i(0) ~ 0.18 for B — > Pq transition is much 
smaller than the typical value of 0.65 ~ 0.70 for the B D transition form factor at q 2 = 0. 



IV. ANALYSIS OF B — ► D**M, P** M DECAYS 
A. Factorization 

In the present work we focus on the Cabibbo-allowed decays B — > P**7r(p), D**D^*\ 
£)**£)(*) anc i _^ D**7i(p), where P** denotes generically ap-wave charmed meson. We will 
study these decays within the framework of generalized factorization in which the hadronic 
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TABLE VI. The dimensionless form factors A and Vb,i,2 at q 2 = m 2 for B — > Z)}^ 2 and 



5 — > D^ 2 transitions calculated in the ISGW2 model. 



Transition A V V\ V 2 

B D\' 2 -0.43 -0.18 0.070 0.49 

B -► £>J /2 0.44 -0.43 -0.15 0.33 

5 S -► D^ 2 -0.49 -0.20 0.085 0.57 

B s -► £>J{ 2 0.50 -0.47 -0.17 0.38 



decay amplitude is expressed in terms of factorizable contributions multiplied by the uni- 
versal (i.e. process independent) effective parameters a« that are renormalization scale and 
scheme independent. Since the aforementioned B decays either proceed through only via 
tree diagrams or are tree dominated, we will thus neglect the small penguin contributions 
and write the weak Hamiltonian in the form 



H, 



cff 



V2{ 



v cb v* 



ai(ob)(sc) + a 2 (sb){cc) + V cb V* d cn(cb)(du) + a 2 (db)(cu) + h.c. (4.1) 



with ((71(72) = 9i7/u(l — 75)92- For hadronic B decays, we shall use a\ = 1.15 and a 2 = 0.26. 

Under the factorization hypothesis, the decays B~ — > D**°D~, B — > D** + Dj and 
5° — > D** + 7r _ receive contributions only from the external IU-emission diagram. As stated 
before, the penguin contributions to the first two decay modes are negligible. 

Apart from a common factor of G F V cb V* d /\/2, the factorizable amplitudes for B~ — > 
L>**°7r- read 



,4(5" 
A(B- 



D*(2308) u vr-) = -aJ.K - m 2 Do )F» v ° 



a 2 f Do (m 2 B -mi)Ft(m 2 Do ), 



D 1 (2427)°7r~) = -2(e* ■ p B )\ ai f n V™ 1 {m^m^ sin 9 + V^ 1 (m^m^ costf 



3/2 



1/2 



+ a 2 



F 1 B7r (m 2 3 / 2 )m 3/2/ 3/2 sintf + F 1 B7r (m 2 V2 )m 1/2/ 1/2 cos6> 



-> D , 1 (2420)°7T-) = -2(e* -p s ) 



3/2 

1 ("4)™ D 3 



3/ 



2 COS — Vr 



BD[ 



1/2 

1 ("4)™ D j 



1/2 



sin# 



+ a 2 



Ff*(m 2 D 3/2)m D 3/2f D 3/2 cos 9 - F^ n (m 2 D i /2 )m D i/2f D i/2 sin6» 



- D 2 *(2460)°7r-) = fc(<) + - m^) + &_(<)< • (4-2) 



Note that except B — > D^tt all other modes receive contributions from color-suppressed 
internal PU-emission. The decay rates are given by 

Pc 



T(B -> L>*7r) 



8nmV y 



D»\ 2 , 



15 



T(B - D*n) 



Pc 



87rm|, l 



m B 



\M(B^D* 2 tt)\ 2 , 



(4.3) 



where A(B — > -D^ 71 ") — £ ]ivPbPb M(B — »■ -D2 71 ") an d is the cm. momentum of the pion. The 
p 2L+1 dependence in the decay rate indicates that only s-, p- and ci-waves are allowed in DqTt, 
D\ix and D^ix systems, respectively. The factorizable decay amplitudes for B~ — > p~ and 
B~ — > D*® p~ are (up to a common factor of GpV c bV* d / \/2) 



A(B~ - D*(2308)V) = 2(e*- PB ) [aif p m p F? D °{m 2 p ) + a 2 f Do m Do A Bp (m 2 Do ) 
A(B~ - L> 2 *(2460)V) = arfpmy^elipB - p D2 )x[ih(m 2 p )e^g au (p p ) p (p p ) a 



+ k(m 2 p )6£5$ + b + {m 2 p )(p p ) a {P P )p9 



(4.4) 



The expression for £? — > D 1 p is more complicated. In the absence of the D^ 2 — D± 2 
mixing, one has 



A(B- - £>i(2420) V) = -iaif p m p [(e* p ■ e* D )(m B + m D AV^ (m 2 ) 



(^■PbWd.-Pb) 



2V2 BDl W) 
m B + m Dl 



- ia 2 f Dl m Dl [(e* • e* Dl )(m B + m^f''^) 

- (g p ■ P B ) (g gl ■ P B ) ~ —7 + ^^uape^e Dl P B PD 1 — 



ra s + Tip 



(4 



Tfl B + m p 

In the presence of the D^ 2 —D^ 2 mixing, it is more convenient to express the decay amplitude 
as 



A[B~ - Dl{e Dl ,p Dl )p-{e p , Pp )] oc e^[S l9uv + S 2 (p B )> B )„ + iS z e, va pp a Di p% (4.6) 

where e 0123 = +1 in our convention, the coefficient S3 corresponds to the p-wave amplitude, 
and Si, S 2 to the mixture of s- and rf-wave amplitudes: 



Si = ai/ p m p (m B + m D 3/2)vf Dl (m 2 ) cos9 — (m B + rripi^V™" 1 (m 2 ) sin 9 
+ a 2 (m B + m p ) 777,^,3/2/^3/2 Af p (m 2 ^) cos 9 - m D i /2 f ' D i/2A Bf '(m 2 ^) sinfl 

1 ..RD^ 2 , ox „ 1 



3/2 



1/2 



S 2 = aif p m p 



Tfi B + m D 3/2 



V 2 1 (m 2 ) sin 



m B + m D i/2 



1 



+ a 2 - 

m B + m p 

S3 = aif p m p 



™> D 3/2j Di 



f n 3/2A^ p (m 2 D y 2 ) cos 9 - m D i/2f D i/2A2 P (m 2 D i /2 ) sin 9 



1 



m B + m D 3/2 



A BD *'\m 2 ) cos 9 



1 



m B + m D i/ 2 



^ D i /2 (m 2 )sin^ 



+ a 2 



m B + ra. 



m D z^2f D z/2V Bp {m 2 D s,2) cos9- m^^^^V^fra^) sinfl]. (4.7) 
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Then the helicity amplitudes H , H + and if_ can be constructed as 

1 



H = 



2m Dl m p 



m 



Di - m p) S l + ^ B Pc S 2 



H± = Si± m B p c S 3 . 



(4.8) 



For B — > Di(2427)°p , the amplitudes 5*1,2,3 ar e the same as Eq. (4.7) except for the 
replacement of cos 9 — > sin# and sinO — > — cos#. The decay rates read (up to the common 
factor of G 2 F \V cb V* d \ 2 /2) 



with 



r(B - 


-Dip) 


T(B - 




T(B - 


-Dip) 



Pc 



D 



iixm 

Pc 



8iTm 2 B 



■\A(B^D* p)/(e*.p B ) 



J p 



(\H \ 2 + \H + \ 2 + \H. 



{ap 7 c + bpl + cp 3 c ), 



(4.9) 



a = Sm B b 2 + 



c = 5m D2 m p k , 

2 rc™2 2 i 2 , o/™2 2 2\i.i j 2n 



6 = 2m| 5 [6m;m| )2 /i 2 + 2(m^ - m^ 2 - m;)£;6+ + k 2 }. 



(4.10) 



B. Results and discussions 

Given the decay constants and form factors discussed in Sec. Ill, we are ready to study 
the B decays into p-wave charmed mesons. The predicted branching ratios are shown in 
Tables VII and VIII. The experimental results are taken from PDG [15] and Belle [14]. For 
B~ — > D*2^ we combine the Belle measurements [14] 

B(B~ -> D* 2 \-)B(Df -> D + 7i-) = (3.4 ± 0.3 ± 0.6 ± 0.4) x 10" 4 , 

b{b~ -> d;\-)b(d;° -> zr+vr-) = (1.8 ± 0.3 ± 0.3 ± 0.2) x io~ 4 , (4.11) 

to arrive at 

B(B~ -> D* 2 °7i-)B(D* 2 ° -> D*+7T-) = (5.5 ± 0.8) x 10~ 4 . (4.12) 

Using i3(L>2 D + n~, D* + n~) = 2/3 followed from the assumption that the Df width is 
saturated by Dtt and _D*7r, we are led to B(B~ — > D^vr - ) = (7.8 ± 1.4) x 10~ 4 as shown in 
Table VII. 

From Table VII we see that except j Di(2427)°tt- the predictions of B(B~ -> D**^-) 
agree with experiment. It is worth mentioning that the ratio 

B(B- -> g|(2460)Vj 

-> £>J(2420)°7r-) 1 J 
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TABLE VII. The predicted branching ratios for B -► D**°(tt- , p~ , D s ) and B~ -► £>**£>(*) 
decays, where -D** denotes a generic p-wave charmed meson. Experimental results are taken from 
PDG [15] and Belle [14]. The numbers in parentheses for D* s0 D and D S \D productions are fits to 
the central values of the experimental data. The axial-vector meson mixing angles are taken to be 
= 17° for D\ — D[ systems and 9 S = 7° for D s \ — D' sl systems. The parameters are taken to 
be ai = 1.15 and a 2 = 0.26. 



Decay 




This Work 


KV [36] 


CM [37] 


KLO [38] 


Expt 




/ > 

D 




DS(2308)°7r- 


7.7 x 10~ 4 


4.2 x 10~ 4 






(9.2 ±2.9) x 10~ 4 


[14] 


R — 




L»i(2427) vr- 


3.6 x 10~ 4 


2.4 x 10~ 4 






(7.5 ± 1.7) x 10- 4 


[14] 


R — 




£>i(2420)°7r- 


1.1 x 1(T 3 


2.1 x 10~ 3 






(1.0 ±0.2) x 10~ 3 
(1.5 ±0.6) x 10~ 3 


[14] 
[15] 


R — 
JD 




D|(2460)°7r- 


6.7 x 10" 4 


7.2 x 10~ 5 


4.1 x 10" 4 


3.5 x 10~ 4 


(7.8 ± 1.4) x 10~ 4 


[14] 


R — 
JD 




£>5(2308)V~ 


1.3 x 1(T 3 












R — 
JD 




L>i(2427)°p- 


1.1 x 1(T 3 












R — 
JD 




L>i(2420)°p- 


2.8 x 10~ 3 








< 1.4 x 10~ 3 [15] 




JD 




L>|(2460)°p- 


1.8 x 1(T 3 




1.1 x 10~ 3 


9.8 x 10" 4 


< 4.7 x 10- 3 [15] 




B~ 




L>j5(2308) l)7 


8.0 x 10" 4 


2.7 x 10~ 3 










B~ 




L>i(2427)° J Dj 


9.6 x 10" 4 


1.4 x 10~ 3 










B~ 




Di(2420)° J Dj 


1.3 x 10~ 3 


5.0 x 10~ 3 










B~ 




D|(2460)°5j 


4.2 x 10~ 4 


1.0 x 10~ 4 


2.7 x 10~ 4 


4.9 x 10" 4 






B~ 




Dj5(2308)°l)*- 


3.5 x 10" 4 












B~ 




Di(2427)°5*- 


6.0 x 10- 4 












B~ 




L>i(2420)° J D*- 


1.6 x 10~ 3 












B~ 




D^(2460)°5*- 


1.1 x 10~ 3 


1.0 x 10" 4 


1.1 x 10~ 3 


1.2 x 10~ 3 






B~ 




D: (2317)- J D 


(9.3 x 10~ 4 ) 









see text 




B~ 




£> a i(2460)-£>° 


(3.1 x 10~ 3 ) 


3.5 x 10~ 3 






see text 




B~ 




D' sl (2536)- D° 


1.3 x 10" 4 


3.4 x 10~ 3 










B~ 




D* s2 (2572)-D° 















B~ 




D* s0 (23l7y D*° 


5.0 x 10" 4 












B~ 




D sl (2460)-D*° 


1.2 x 10~ 2 












B~ 




D' sl (2536)- D*° 


5.3 x 10~ 4 












B~ 




D* s2 (2572)-D*° 
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TABLE VIII. Same as Table VII except for neutral B and B s mesons. 



Decay 


This Work 


KV [36] 


CM [37] 


KLO [38] 


Expt 




S u -► D5(2308)+7T- 


2.6 x 1(T 4 


4.1 x 10" 4 










B° -► L>i(2427) + vr" 


6.8 x 1(T 4 


1.2 x 10" 4 










B Di(2420)+7r- 


1.0 x 10" 3 


2.4 x 10" 3 










S° D|(2460)+7r- 


6.1 x 1(T 4 


7.1 x 10" 5 


4.1 x 10" 4 


3.3 x 10" 4 


< 2.2 x 1(T 3 


[15] 


S u -► L>j$(2308)+p" 


6.4 x 10" 4 












5° -► L>i(2427)+p- 


1.6 x 10" 3 












5° -► Di(2420)+p- 


2.6 x 1(T 3 












B° L>?j(2460) + p" 


1.7 x 10" 3 




1.1 x 10" 3 


9.2 x 1(T 4 


< 4.9 x 10 -3 


[15] 


S U -► D5(2308) + L»7 


7.3 x 10" 4 


2.6 x 10" 3 










5° -► L>i(2427)+£> s " 


8.8 x 1(T 4 


1.3 x 1(T 3 











B° 
B° 



Di(2420)+ D; 
D* 2 (2460)+D- 



1.2 x 10" 3 
3.8 x 10" 4 



4.9 x 10" 3 

1.0 x 10" 4 2.7 x 10" 4 4.6 x 10" 4 



W 
B° 
B° 
B° 



Z^(2308)+L> 
L>i(2427)+£>* 
D[(2420)+ D* s 
L>?j(2460) + Z)* 



*— 

s 



3.2 x 
5.5 x 
1.5 x 
1.0 x 



10" 4 
10" 4 
10" 3 
10" 3 



1.0 x 10" 4 1.1 x 10" 3 1.1 x 10" 



W 

B° 
B° 
B° 



D} 
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is measured to be 0.77±0.15 by Belle [14] and 1.8±0.8 by CLEO [18]. The early prediction by 
Neubert [39] yields a value of 0.35. Our prediction of R = 0.61 is in accordance with the data. 
However, the predicted rate for D 1 (2427)°ir~ is too small by a factor of 2. This is ascribed 
to a destructive interference between color-allowed and color-suppressed tree amplitudes 

because the form factors V 1 and V 1 have signs opposite to that of F± as required 
by heavy quark symmetry (see Eq. (4.2) and Table VI). In contrast, the production of 
_Di(2427) + 7r~ is larger than .D 1 (2427)°7r~ by a factor of about 2 because the former does 
not receive the destructive contribution from internal iy-emission. Hence, a measurement 
of the ratio _Di(2427)°7r"/.Di(2427) + 7r~ can be used to test the relative signs of various 
form factors as implied by heavy quark symmetry. Note that for the D 1 ^ 2 — D^ 2 mixing 
angle we use 9 = 17° [see Eq. (2.7)]. If a negative value of —17° is employed, the decay 
B~ — > /} 1 (2427)°7r~ will be severely suppressed with a branching ratio of order 6 x 10~ 6 . This 

1/2 3/2 

means that the D] 1 - D\' mixing angle is preferred to be positive. In their study Katoch 
and Verma [36] obtained a small branching ratio for B~ — > Dq(2308) 7t~ as they assumed 
a vanishing decay constant for Dq. As stressed before, this decay constant is comparable to 
fo because of charm and light quark mass imbalance and a rough estimate yields fp* ~ 160 
MeV [cf. Eq. (3.5)]. Consequently, the contribution from internal iy-emission will account 
for the aforementioned discrepancy between theory and experiment. Moreover, the ratio 
Dq + ii~ /Dq°it~ is predicted to be 0.34 instead of unity because of the absence of the color- 
suppressed tree contribution to the former. 

At a first glance, it appears that the prediction B(B~ -> L>i(2420)°p~) = 2.8 x 10~ 3 
already exceeds the experimental limit 1.4 x 10~ 3 [15]. However, it should be noticed that the 
D[°p~ rate is about three times larger than that of D'®tt~ as expected from the factorization 
approach and from the ratio f p /f n ~ 1.6 (see Table IX below). Hence, it appears that the 
present limit on D[p~ is not consistent with the observed rate of D\-k~ . Of course, it is 
crucial to measure B — > D' 1 (2A20)p in order to clarify the issue. 

Apart from the external H^-emission diagram, the D**n(p) productions in neutral B 
decays also receive H^-exchange contributions which are neglected in the present work. This 
will constitute a main theoretical uncertainty for B° — > D** + n~ (p~). 

For B — > DD** decays, their factorizable amplitudes are rather simple 

A(B^DDT) = ^^^(STI^IOXDK^IS). (4.14) 

The decay rates are thus governed by the decay constants. Since the tensor meson cannot 
be produced from the V — A current, the B decay into D* 2 D^ is prohibited under the 
factorization hypothesis. For D* Q and D sl mesons, it is expected that fjj* ~ /d s i(2460) ^ 
Id' (2536) as fn 1 / 2 = Id* and / n i/2 = in the heavy quark limit [see Eq. (3.7)]. When the 
charmed quark mass is finite, the Dl^—D^J 2 mixing effect has to be taken into account via 

/d.i = /nVa COS 9 S + / 3/2 SU1 a , 

f D > = -/ D i/2 sin6» s + / 3/2 cos6» s , (4.15) 
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where use of Eqs. (2.8) and (3.11) has been made. Consequently, it is anticipated that 
DD* s0 > DD s i(2A60) 3> DD sl (253Q) after taking into account phase space corrections. How- 
ever, the recent data on DD** production indicate a surprise. 
The results of Belle measurements read [40] 

B[B -> D J D s * (2317)]i3[ J D: (2317) -> D s n } = (8.5i?;J ± 2.6) x 10" 4 , 
B[B -> D£> s i(2460)]B[Z> s i(2460) -> £>>°] = (17.8l£;jj ± 5.3) x 10" 4 , 
B[£ -> DD s i(2460)]B[D fl i(2460) -> ^7] = (6.7±^;1 ± 2.0) x 10" 4 . (4.16) 

Since -D* (2317) is dominated by its hadronic decay to D s -k ^ the branching ratio of B — > 
-D.D* (2317) is of order 0.9 x 10~ 3 . This means that the production rate of D°D* in £? decays 
is smaller than B(B + — > £) D s ) = (1.3 ± 0.4)% [15] by one order of magnitude. Since this 
decay proceeds only via external l/F-emission, it can be used to determine the decay constant 
of D* s0 . It is found that f D * ~ 60 MeV M which seems too small, recalling that fo* ~ 160 
MeV is needed to account for the production of -Do 71 "- To estimate the branching ratios 
of D*7T° and in .D sl (2460) decay, we need some experimental and theoretical inputs. 
Besides the values 0.38 ± 0.12 and 0.63 ± 0.21 measured by Belle for the ratio -D s 7/-D*7r° 
[see Eq. (2.20)], the ratios D*j/D*Tt° and D s n + 7i~ / D*n are found to be less than 0.16 and 
0.08, respectively, by CLEO [13]. Theoretically, the Ml transition D s i — > -D* 7 turns out 
to be quite small [41]. Assuming that the -D s i(2460) decay is saturated by D* s tt°, D s 7, D* s ^ 
and D s -K-n, we are led to 

0.56 < £(£> sl (2460) -> D* s n°) < 0.75, (4.17) 

where the averaged value 0.44 ± 0.10 for the ratio D s / ~f/D*7i° has been employed. This in 
turn implies B[B — > DD s i(24Q0)] = (1.5 ~ 4.4) x 10~ 3 . As a result, the decay constant of 
£> s i(2460) is found to be 

f Dsl (2460) ~ 100 - 175 MeV. (4.18) 

This is in sharp contrast to the heavy quark symmetry expectation that the decay constants 
of D* s0 and D sl become identical in the heavy quark limit. Therefore, the large disparity 
between fu* and fo sl is a surprise. It is not clear to us what is the underlying reason for 
the discrepancy between theory and experiment. 

Although the B decay into D* s2 D^ is prohibited under the factorization hypothesis, 
nevertheless it can be induced via final-state interactions (FSIs) and/or nonfactorizable con- 
tributions. For example, it can be generated via the color-allowed decay B~ — > D S D^ 



ttThe upper limit on the ratio T(D* s0 -► D* s j)/T(D* s0 -► D s 7r°) < 0.059 was set recently by CLEO 
[13]. 

^Interestingly, this happens to be the original estimate of fo* made in [4] for D* in the four-quark 
state. 
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followed by the rescattering process D S D^ — > D s2 D°. Since the nonfactorizable term is of 
order c-ijN c with the Wilson coefficient c^ipi^) ps —0.20, it is likely to be suppressed relative 
to the FSIs. Hence, an observation of B — > D s2 D^ could imply the importance of final-state 
rescattering effects. 

Since heavy quark symmetry implies /d' 1 (2536) ^ /r> s i(2460), it is important to measure 
the B decay into D sl (2536)D^ to see if it is suppressed relative to D s i(24:60)D^ to test 
heavy quark symmetry. In order to compute B — > D s D* we have fixed the decay constants 
of .D* , D l J 2 and D 3 S { 2 as shown in Eq. (3.11) so that the resulting production rates of D* s0 D° 
and D s iD° are consistent with experiment (see Tables VII and VIII). We see that the decay 
B — > D s i(2A60)D*° has the largest branching ratio of order 10 -2 in two-body hadronic B 
decays involving a p-wave charmed meson in the final state. It is essential to test all these 
anticipations in the near future. A recent CLEO measurement yields [42] 

B[B~ -> (D s + £>*)(£>? + D[° + Df)} = (2.73 ± 0.78 ± 0.48 ± 0.68)%. (4.19) 

Our prediction of 6 x 10~ 3 is slightly small. 

It is interesting to consider the ratios T(B -> D**V)/T(B -> D**P) and T(B -> 
D**y)/r(-B — > D**P). The calculated results are shown in Table IX. Several remarks 
are in order: (i) The ratios D**° p~ / D**°n- and D** + p~ / D** + 7r~ for D** = D* and D x are 
not the same as the former receives an additional color- suppressed internal VT^-emission con- 
tribution, (ii) Whether the ratio D**V/D**P is greater than unity or not depends essentially 
on the ratio of the decay constants. For example, DqD*/DqD s = 0.43 if the decay constants 
of D* and D s are similar, while Dq + p~ / Dq + h~ = 2.5 for f p /f w = 1.6 . It should be stressed 
that the proximity of the ratio D** + p~ / D** + n~ to 2.5 has less to do with the three degrees of 
freedom of p, rather it is mainly related to the decay constant ratio of f p /f n - (hi) The ratio 
of D**p~ /D**ir~ in B s decay is the same as D** + p~ / D** + n~ in B decay as they proceed via 
external iy-emission. 

TABLE IX. The ratios F(B -► D**V)/T{B -► D**P) and T(B D**V)/T(B D**P). The 
last column is for the ratio T(B S D**p-)/r(fi s -»• £>**tt-). 
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D**°p- 

£)**U 7r - 
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0.43 
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0.63 
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2.5 




2.6 


2.6 


1.2 


4.0 


2.6 




2.7 


2.7 


2.7 




2.7 



Because the scalar resonances Dq and Di have widths of order 300 MeV, we have checked 
the finite width effects on their production in B decays and found that the conventional nar- 
row width approximation is accurate enough to describe the production of broad resonances 
owing to the large energy released in hadronic two-body decays of B mesons. 



22 



C. Comparison with other works 



The decays B — > D**(tt,D s ) and B — > D**D have been studied previously by Katoch 
and Verma (KV) [36]. Lopez Castro and Muhoz (CM) [37], Kim, Lee and Oh (KLO) [38], 
and Lee [43] also have a similar study with focus on the tensor charmed meson production. 
We shall comment on their works separately. 

In the paper of KLO, the B — > D 2 form factors are evaluated using the ISGW2 model. 
However, the predicted rates for D 2 ir and D 2 p by KLO are smaller than ours by a factor 
2 (see Tables VII and VIII); that is, their B — > D 2 form factor differs from ours roughly 
by a factor of In contrast, the results for D 2 D^f> are similar owing to the fact that 
KLO employ fu s = 280 MeV and fu* = 270 MeV which are larger than ours [see Eq. 
(3.11)]. Instead of employing the ISGW2 model for B — > D2 form factors, Lee applied the 
QCD sum rule to study the universal function t 3 / 2 (c<j) to the leading order [43]. He found 
B(B~ — > D 2 °tt~) = 12.9 x 10~ 4 for a\ = 1.15, which is two times larger than ours. 

The work of KV and CM is based on the original version of the ISGW model. However, 
as stressed by KLO [38], the exponentially decreasing behavior of the form factors in the 
ISGW model is not realistic and justified. This has been improved in the ISGW2 model 
which provides a more realistic description of the form-factor behavior at large (q^ — q 2 ). 
The values of form factors at small q 2 in the ISGW2 model can be a few times larger than 
that obtained in the ISGW model as the maximum momentum transfer q^ in B decays is 
large. However, the expected form-factor suppression does not appear in the calculations of 
KV as they calculated the form factors at q^ and then employed them even at low q 2 . In 
contrast, CM did compute the form factors at proper q 2 . The fact that CM and KLO have 
similar results for B -> D* 2 it(p) and B -> D*D^ (see Tables VII and VIII) is surprising as 
the B — > D 2 transition is evaluated in two different versions of the ISGW model. To check 
this, we find that i](q 2 ) = k(q 2 ) + b+(q 2 )(m 2 B — m 2 D2 ) + b_(q 2 )q 2 = 0.286 and 0.386 for q 2 = ml 
in the ISGW and ISGW2 models, respectively, which in turn imply the respective branching 
ratios, 6.7 x 10~ 4 and 3.8 x 10~ 4 for B~ — > D 2 °7r~ . Therefore, our result obtained in the 
ISGW model is consistent with that of CM and the estimate of f](q 2 ) within the ISGW2 
model by KLO is likely too small, as noted in passing. 

The axial- vector charmed meson production considered by KV is for Di( l Pi) and Di( 3 Pi) 
rather than D^ 2 and D^ 2 . Therefore, the expressions of the decay amplitudes involving D\ 
or D s i by KV should be modified by taking into account a proper wave function combination, 
Eq. (2.1). For the decay constants, KV assumed that f D * = and fo^p^ =0. As a 
consequence, B~ — > Dq°ti~ is too small compared to experiment and the decay into D s0 D 
is not allowed. This is not consistent with the heavy quark symmetry relation f n i/2 = fo* ■ 
Finally, the predicted rate of B~ — > D 2 tt~ by KV is too small by one order of magnitude 
compared to experiment. This is ascribed to a missing factor of (ms/mr) 2 in their calculation 
of decay rates. 



23 



V. CONCLUSIONS 



The hadronic decays of B mesons to a p-wave charmed meson in the final state are studied. 
Specifically we focus on the Cabibbo-allowed decays B — > D**7i(p), D**D^*\ D**D^ and 
B s — > D**n(p). The main conclusions are as follows: 

• We apply heavy quark effective theory in which heavy quark symmetry and chiral 
symmetry are unified to study the strong decays of p-wave charmed mesons and de- 
termine the mag nitude of the D\ l2 -Dl 12 mixing angle. In contrast, the present upper 
limits on the widths of -D sl (2460) and .0^(2536) do not provide any constraints on 
the D l J 2 — D z J 2 mixing angle 9 S . Therefore, we appeal to the quark potential model to 
extract 9 S . 

• Various form factors for B — > D** transitions and their q 2 dependence are studied 
using the improved version of the Isgur-Scora-Grinstein-Wise quark model. Heavy 
quark symmetry constraints are respected in this model calculation. 

• The predicted branching ratios for B~ — > D**ir~ agree with experiment except 
Di(2427)°7r~ . The D\^ 2 —D\^ 2 mixing angle is preferred to be positive in order to avoid 
a severe suppression on the production of Di(2427)°7r _ . The decay B~ — > _D / 1 (2420)°p~ 
is predicted at the level of 3 x 10~ 3 . Although it exceeds the present experimental limit 
of 1.4 x 10 -3 , it leads to the ratio D / 1 (2420)p~/D / 1 (2420)tt~ « 2.6 as expected from the 
factorization approach and from the ratio f p /f n ~ 1.6 . Therefore, it is crucial to have 
a measurement of this mode to test the factorization hypothesis. 

• The predicted rate for B~ — > /}i(2427)°7r~ is too small by a factor of 2 owing to a 
destructive interference between color-allowed and color-suppressed tree amplitudes as 

1/2 3/2 

the relevant form factors for B — > D x and B — > D\ transitions are negative. It is 
crucial to measure the production of .D 1 (2427) + 7r~ to see if it is larger than .D 1 (2427)°7r~ 
by a factor of about 2 because the former does not receive the internal lU-emission 
contribution. 

• Under the factorization hypothesis, the production of D s2 D^ in B decays is prohibited 
as the tensor meson cannot be produced from the V — A current. Nevertheless, the 
decays B — > D* s2 D^ can be induced via final-state interactions and/or the nonfactoriz- 
able contributions. Since the latter is suppressed by the order of C2/N c , an observation 
of B — > D* s2 D^ could imply the importance of final-state rescattering effects. 

• For B — > D s D decays, it is expected that D* s0 D > D s iD as the decay constants of 
the multiplet {D* , D sl ) become identical in the heavy quark limit. The experimental 
observation indicates that D* s0 D/D sl D w (0.2-0.6) and that f Dsl m (100- 175) MeV 
and f D * ~ 60 MeV. The less abundant production of D* s0 D than D sl D and the large 
disparity between f Dsl and fo* are surprising. The reason for the discrepancy between 
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theory and experiment remains a puzzle. In the meantime, it is also important to 
measure the decay to D sl (2536)D^ to see if it is suppressed relative to D sl (2A60)D^ 
to test the heavy quark symmetry relation /d s i(2536) < /r> s i(246o)- 
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Note Added: 

After this work was completed, we noticed the appearance of the related works on B — > DD** 
decays by C.H. Chen and H.n. Li, hep-ph/0307075; A. Datta and P.J. O'Donnell, hep- 
ph/0307106; M. Suzuki, hep-ph/0307118. 
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